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Abstract 



The dispersion of one hole in an extended t-J model with additional hopping terms 
to second and third nearest neighbours and a frustration term in the exchange part has 
been investigated. Two methods, a Green's function projection technique describing 
a magnetic polaron of minimal size and the exact diagonalization of a 4 * 4 lattice, 
have been applied, showing reasonable agreement among each other. Using additional 
hopping integrals which are characteristic for the CuC>2 plane in cuprates we find in 
the nonfrustrated case an isotropic minimum of the dispersion at the point (n/2, n/2) 
in fc-space in good coincidence with recent angle-resolved photoemission results for the 
insulating compound S^CuO^C^. Including frustration or finite temperature which 
shall simulate the effect of doping, the dispersion is drastically changed such that a 
flat region and an extended saddle point may be observed between (w/2,0) and (tt,0) 
in agreement with experimental results for the optimally doped cuprates. 



PACS numbers: 74.25. Jb, 74.72.-h,75.10.-b 



1 Introduction 



It is well known that many properties of the high-T c superconductors (HTSC's) are deter- 
mined by the spectrum of holes in the doped two-dimensional (2D) Cu02 planes. Recent 
experiments by angle resolved photoemission (ARPES) reported a small band width pro- 
viding evidence for the existence of strong correlations in cuprates [1-5]. Moreover they 
indicate the existence of a flat dispersion region close to the Fermi surface for the optimally 
doped cuprates Bi2212, Bi2201, Y123 and Y124 (for the notation see g). The flat region 
of the quasiparticle band is located close to the point (it, 0) in the irreducible part of the 
Brillouin zone (where the lattice constant has to be set to unity a = 1) and it has the form 
of an extended saddle point between (vr, 0) and (tt/2, 0). The existence of the flat band 
region leads to the appearence of a strong Van-Hove (VH) singularity which is close to the 
Fermi level. This circumstance is often used by different theories for the explanation of the 
high superconducting transition temperature [3, 6-9]. Then, the existence of the optimal 
value of doping is naturally connected with the movement of the chemical potential over 
the position of the VH singularity when the hole density is increased. So, the investigation 
of the hole spectrum in the Cu02 plane is of great importance for a microscopic theory of 
HTSC. 

Important new results were obtained recently by ARPES experiments for the insu- 



lating, antiferromagnetic compound Sr2Cu02Cl2 [1C]. These results correspond to the 
situation of one hole in the quantum antiferromagnet and they serve as a test for different 
theoretical approaches. It was found p3] that the bottom of the hole band is close to the 
point (vr/2, vr/2) and the band width is of order 0.3 eV. It is essential, that the effective 
mass was shown to be isotropic at the bottom of the band. 

One generally believes that the t-J model describes qualitatively the hole spectrum 
of the Cu02 plane in HTSC. There are numerous studies of the hole spectrum in the 
framework of the t-J model based on the exact diagonalization of small clusters ]ll], 12], 
the self-consistent Born approximation (l^, |l4|] or a "string" ansatz for the hole wave 
function ipHft . One can either start from the two-sublattice Neel-type state [jD^, 14, |l5[ or 
a spin rotational invariant spin liquid state juj and obtains qualitatively the same result: 
the hole motion occurs mainly on one sublattice, i.e. the dispersion relation is dominated 
by an effective hopping to next nearest neighbours with the minimum of the dispersion 
at (7r/2, 7r/2). One finds a flat region of the dispersion near to the point (tt,0). That was 
used for an attempt to interpret the experimental data of the optimally doped compounds 



[17], where a rigid-band picture was assumed (see also Iflql). However, the anisotropic 
minimum of the pure t-J model is in contrast to the experimental result in Sr2Cu02Cl2 
which shows an isotropic band bottom and a large energy difference between (tt/2, 7r/2) 
and (-7T, 0). 

Recently it was shown [19] that the qualitative agreement with the experimental results 
|l0| may be improved in the framework of the t\-t2-J model, which takes into account the 
next-nearest-neighbour hopping ti of the hole. The inclusion of a t2-term substantially 
improves the correspondence of the calculated spectrum near the band bottom with the 
experimentally observed isotropic minimum of e(k) at the point (tt/2,tt/2). On the other 
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hand, the hole spectrum of the tx-t^-J model does not reproduce the flat band region near 
the point (vr,0). So, we have the following theoretical problem: the one-hole approach in 
the pure t-J model can reproduce the flat band region of the optimally doped cuprates [1- 
5] but not the isotropic band bottom of the insulating compound [10]. And vice versa, one 



hole in the tx-t^-J model leads to an isotropic band bottom but not to a flat band region 
near (tt, 0). One possible solution could be that we have to choose different microscopic 
models for the different compounds. Here, we will investigate the other possibility, whether 
it is possible that the hole band is changed by doping so that a flat band region between 
(tt, 0) and (tt/2,0) arises. 

To simulate the effect of doping in a study of the one-hole motion there were recently 
compared two possibilities, namely the influence of frustration and temperature p0|. The 
inclusion of frustration leads to a t-Jx-Ji model, where J\ and Ji denote the antiferromag- 
netic exchange interaction between nearest and next nearest neighbours. The influence 
of frustration was investigated by two different methods, namely by a variational ansatz 
where the spin-spin correlations in the frustrated Heisenberg model had been calculated 
by a spin rotational invariant procedure pi}] and by the exact diagonalization of a 4 * 4 
lattice. It was found that the frustration shifts the minimum of the hole dispersion from 
(tt/2, 7t/2) to the point (tt, 0). Both methods had a quite reasonable agreement where the 
variational method showed the effect in a more pronounced way. Flat dispersion regions 
were found in the nonfrustrated and frustrated cases. Without frustration, the flat region 
occurs around (7r/2,7r/2), whereas for frustration J2/J1 ~ 0.4 it occures between (vr,0) 
and (tt, 7r/2). Such a flat region is similar to the experiment [1-5], but there it appears 
between (vr,0) and (n/2,0). 

The mentioned above investigations demonstrate that the extension of the t-J model 
by the inclusion of additional hopping terms, frustration or by taking into account a finite 
temperature may strongly modify the initial features of the hole spectrum. But at the 
same time, none of these extensions alone can lead to an adequate description of the 
experimental picture. 

In the present paper we will investigate the hole spectrum in the framework of the 
t\-t2-t^-J\-J2 model, e.g. we will take into account the hole hoppings between first (t\), 
second (£2) and third (£3) nearest neighbours and the frustration in the spin subsystem 
simultaneously. The additional hoppings £2 and £3 naturally appear in the Hamiltonian if 
one reduces the well known three band Hubbard model of the CuC>2 plane to an extended 



t-J model p2|, 23, M. Throughout the present paper we will fix £2 and £3 to such values 
which can be derived from the characteristic parameters of the three band Hubbard model 
for the Cu02 plane which were given by Hybertsen et al. p5[]. The frustration will be 
understood to simulate the doping as proposed, for instance, in Ref. |p6[] . The effect of a 
finite temperature will also be discussed. 

Two methods will be used to investigate the problem. The projection Mori-Zwanzig 
method [^?J for two-time retarded Green's functions will be employed to treat the hole 
excitation as a magnetic polaron of minimal size. Simultaneously we will use the exact 
diagonalization of a 4*4 lattice. Within the projection method the dispersion is determined 
by static spin-spin correlation functions. These are calculated here in a spin rotational 
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invariant method [21] (see also the zero temperature version |28|| ). We describe the state of 
the magnetic subsystem as a spin liquid state in contrast to the widely used two-sublattice 
Neel-type state due to several reasons. At first, our choice gives the possibility to avoid the 
degeneracy of the hole-spectrum between the points (0, 0) and (ir, it) which always occurs in 
the Neel-type state. Second, the spin liquid state does not contradict the Mermin- Wagner 



theorem [29] which forbids magnetic order in two dimensions for any finite temperature. 
So it gives the possibility to investigate the one-hole motion at non-zero temperatures. 
Of course, the Mori-Zwanzig projection method is an approximate one. Therefore, it is 
reasonable to investigate the problem by the exact diagonalization method and to compare 
the results of two independent procedures. 

The paper will be organized as follows. At first, in Sec. 2 we will present the spin 
polaron approach in the framework of the projection method. Then we show the results 
of both methods for zero temperature and increasing frustration (Sec. 3) and discuss the 
influence of temperature. Finally we present our conclusions. 



2 Small spin-polaron approach for the hole spectrum 

We consider an extended t-J model on a square lattice with hoppings between first, second 
and third nearest neighbours and with a frustration term in the exchange interaction. The 
Hamiltonian of the model is given by 

H = H t + Hj = — ^ tl X(°X^i + - '^Jm Si ■ Si+m , (1) 
ial im 



where 



ti > if I is a nearest neighbour vector 

ti if I is a second nearest neighbour vector 

t^ if I is a third nearest neighbour vector 

else 



and 



J\ > if m is a nearest neighbour vector 
Jm= \ J2 > if m is a next nearest neighbour vector . (3) 
else 

The Hamiltonian is expressed in terms of the Hubbard projection operators which exclude 
the double occupancy at site i: acting on the vacuum state X[° creates the electron 
(annihilates the hole) at site i with spin S = 1/2 and spin projection a/2 (a = ±1 is a 
spin index). The following operator equations are fulfilled 

Xf = 4(1 - tw) , = n la (l - tw) , Af- CT = cla^c , 

nia = c\ a c la , SI = Xr a , S? = \ £ CT aXr , (4) 

xf° + xr = 1 , X^X^ = X^5 X2X3 , 

where A n = or a in the last equation. The hopping parameters t\ in (Q) are not free 
but have to be determined by a reduction procedure from the three band Hubbard model 
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of the C11O2 plane. We use here an analytical reduction procedure, the so-called "cell 



perturbation method" [22, 23, E4[]. For the parameters of the three band Hubbard model 
we choose the values which were found by Hybertsen et al. |25| for La2Cu04 by a 
constrained density-functional calculation, namely: t p d = 1.3 eV, t pp = 0.65 eV,e p — £d = 
3.6 eV, Ud = 10.5 eV, U p = 4 eV and U p d = 1.2 eV |Pq] . Band structure calculations for 



Sr2Cu02Cl2 [ 30 1 show only few differences to the La-compound. Therefore, we expect 
that the parameter values of Hybertsen et al. are also responsible for the Sr-compound. 
From the reduction procedure (details will be published in Ref . p3] ) we found the following 
hopping parameters: t\ = 498 meV, ti = —41 meV and t% = 77 meV. They are within the 
limits which have been found in [^] and also near to the results of a numerical mapping 
in ||25}|. But in difference to Hybertsen et al. we found it necessary to include the t3-term 
also. In the following we will choose t\ as the unit of energy, i.e. t% = —0.08, £3 = 0.15. So 
you see that the additional hopping terms are rather small, but nevertheless essential as 
will be shown below. The reduction procedure gives also the possibility to calculate the 
exchange energy J\. That is however less important for our purpose, since we will show 
that the results for different values of J\ can be roughly transformed into each other by 
scaling the bandwidth with J\. 

Let us discuss the one-hole excitations on the background of the half-filled rotationally 
invariant singlet states |^) of the pure spin system Hj. We are interested in the energet- 
ically lowest branch of these excitations and it is known that the bottom of its band is 
represented by the states with total spin Stot = 1/2- We will treat the problem within the 
framework of the spin polaron of small radius. That means that we restrict ourselves to 
spin excitations in the immediate neighbourhood of the hole. It may be seen that the full 
basis of such excitations is given by the following 5 basis operators (j>f* with a = 0, 1, ... 4: 

4 = X? , tf = -£X?I a X?° (a = l,...,4), (5) 

s 

where we will use the following notation hereafter: the small latin letters a and b denote 
either a number between and 4 or the corresponding lattice vector 

0^(0,0) 1^(1,0) 3 ~(-l,0) , . 

2^(0,1) 4^(0,-1) {) 

in a synonymous way. Any vector (f>^ \^>) corresponds to a one hole state with Stot = 1/2 
and Sf ot = —a/2, where S tot and St ot are the spin and its projection of the total system. 

To obtain the one-hole spectrum e(k) we use the two-time retarted matrix Green's 
function G ab (t,k) for the Fourier transformation of the operators 

G a \t,k) = mty^k)) = ^ f (o)» , (7) 

*8 = -^E e<iy #. ( § ) 



N , 



where {...,...} denotes the anticommutator and where we have used Zubarev's notation 
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The equation of motion for the Green's function (|?]) after going over from time t to 
frequency u by Fourier transformation {{4>i; (fy })<j has the form 

u{{<t>l\<t> h k))u = S ab (k) + M;^))^ (9) 
S ab (k) = m,^}), <fl = {<t>lH], (10) 

with the comutator [...,...]. To restrict ourselves to the above chosen set of 5 operators 
{4>%} (|D we use the standard Mori-Zwanzig projection method p7| for the operators ip%: 

<p% ~ ^O" 6 ^)^; m=mk)S-\k); (11) 

H a \k) = ({[^,fl],A = (f|0^V)-rt), (12) 

where Eq denotes the ground state energy of the spin system Hj without any hole. Note 
that the expression with the anticommutator in (|l^) can be simplified considerably for the 



one-hole problem. Using (12) the eqn. (||) takes the following matrix form 

(ojE - H(k)£r 1 (k)) G = S{k) , (13) 



where E_ is a unit matrix. The quasiparticle spectrum e{k) is determined by the poles of 
the Green's function G from the condition 

det\e(k)S{k)-H(k)\ = . (14) 

The matrix elements of H_ and S_ may be calculated in a straightforward way and the 
explicit expressions are given in the Appendix. These matrix elements are expressed 
through spin-correlation functions of the spin system for an undoped lattice which is 
described by the frustrated 5 = 1/2 Heisenberg model. As mentioned in the Introduction 
we treat this model in the framework of the spherical symmetric approach which gives 
values of pair spin-spin correlation functions at any value of frustration parameter J2 / J\ 
and temperature. Let us note that the Ritz variational principle with the 5 basis vectors 
(f)f* 1^) leads exactly to the same condition ( |i~4|) and the same matrices H_ and S_ |f[6| , |2"of |. 
So, at least for the case of the one-hole motion, the Mori-Zwanzig projection technique is 
identical to the Ritz variational principle. 



3 Numerical results and discussion 



As mentioned above we calculate the spectrum e(k) using the hopping parameters t\, £2, 
£3 which were obtained for La-cuprates and take t\ as the unit of energy. So we have the 
values t2 = —0.08 and £3 = 0.15. In the calculations of the magnetic polaron of minimal 
size we choose J\ = 1 since our relevant set of operators is too limited to describe the 
small Ji case. But we will see that one can reach the physically relevant region by a simple 
rescaling. 

Simultaneously with the Green's function method the Lanczos exact diagonalization 
scheme is used as a complementary method to calculate the low-lying eigenstates for a 
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square lattice of 4 * 4 sites (with periodic boundary conditions). These low-lying states 
are classified by its momentum k and the total spin. Here we will concentrate on the 
band with total spin 1/2. There are only few exceptions where the lowest states have spin 
3/2 and that happens only for higher lying levels with momentum (7r,7r) or (0,0). The 
Mori-Zwanzig projection method gives the energy difference between the lowest state with 
one hole for total spin 1/2 and fixed momentum and the ground state without any hole. 
The same energy difference is calculated for the 4*4 lattice. 

To clarify the importance of the additional hopping terms and frustration, in Fig. 1 we 
represent the results for J\ = t\ = 1 without frustration (Fig. la) and with J2 = 0.4 (Fig. 
lb). For convenience we have chosen a finite temperature T = 0.2 to calculate the static 



spin-spin correlation functions [21, 28], but we have checked that the differences to the zero 
temperature case can be neglected. One observes a quite reasonable agreement between the 
Green's function method and the exact diagonalization data for these parameter values. As 
it is seen in Fig. la, in the absence of frustration, the spectrum demonstrates an isotropic 
band bottom close to the point (jv/2, 7r/2) in accordance to the experimental dispersion of 
one hole |10| . The same result was obtained in Ref. |l9| for the t\-t2-J model where the 
following parameter values were taken: £2 = — 0.35ii, J = 0.3t\. Let us mention that the 
isotropization of the band bottom in comparison with the t-J model may be obtained in a 
rather wide range of parameters £2 and £3. For example, also the values £2 = — 0.15ii and 
£3 = O.lti give a dispersion which is similar to that one shown in Fig. la. It is important 
that the spectrum in Fig. la does not demonstrate a flat band region near the point (tt, 0). 

Before considering the spectrum of Fig. lb let us remind that we suppose that the 



frustration simulates doping [26|. Of course there is no full equivalence between doping 
and frustration. For example, the doped t-J model and the frustrated J1-J2 model give 
different results for the dynamical spin-spin structure factor and for the spectrum of Raman 



scattering [32|. Nevertheless, it is well known that both doping and frustration lead to a 
decrease of the magnetic correlation length. Furthermore, numerical calculations on finite 
lattices indicate the equivalence of the mentioned models if we are interested in the static 



spin-spin correlation functions [33|. Note that this is especially relevant in our present 
Green's function method where the spectrum e(k) is determined by the static spin-spin 
correlation functions of the spin subsytem. 

Fig. lb represents the spectrum for the same energetical parameters as Fig. la, except 
the inclusion of a frustration term J2 = 0.4 Ji . The comparison of Figs, la and b indicates 
that the frustration dramatically changes the spectrum in the vicinity of the point (tt, 0). 
It leads to the appearence of a flat band region close to that point. Moreover, this flat 
band region has the form of an extended saddle point which stretches in the direction 
(n/2, 0)-(-7r, 0). Such a band structure corresponds to the ARPES results for optimally 
doped cuprates [1-5]. 

So, we can describe the experimental results for insulating (no frustration, Fig. la) 
and optimally doped compounds (Fig. lb, the case with the frustration) by the same set 
of hopping parameters. It is important that at the same time our approach demonstrates 
the non-rigid band behavior in the framework of a simple and natural mechanism: the 
doping leads to the frustration in the spin subsystem and to the variation of spin-spin 
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correlation functions, and correspondingly the alteration of these functions results in the 
non-rigid behaviour of the spectrum. Or, in other words: the doping changes the state 
of the magnetic background which strongly determines the form of the spectrum. Let us 
mind that this mechanism doesn't take into account the direct interaction between holes 
which will be of course important in the regime of heavy doping. 

In Figs. 2a and b we compare the exact diagonalization and the Green's function 
projection method for several momenta in dependence on the frustration Jij J\. The 
non-rigid band behaviour is clearly seen. Both figures demonstrate the following common 
features. First, without frustration, the order of the levels coincides. Then, the change 
of the energy due to frustration is quite similar for the lowest four levels. The energy 
difference between the levels (tt/2, it/2) and (vr,0) decreases as the frustration increases, 
i.e. the level (tt, 0) goes down with the respect to the band bottom (in the case t<2 = t% = 
there is even an inversion of the levels and the band bottom occures close to (n, 0) pC| ). 
One finds a crossing for the levels (n/2,0) and (ir,ir/2) at J2/J1 = 0.17 and 0.27 for 
the projection and exact diagonalization methods, respectively. In addition, the energy 
difference e(ir, 0)—e(ir/2, 0) decreases in both methods, and both levels cross in the Green's 
function projection technique. That corresponds to the appearance of the flat band region 
and the extended saddle point in the vicinity of the (tt, 0)-point. Of course, there are also 
remarkable deviations of the Green's function method, especially for the level (it, it). But 
fortunately, that concerns only a high lying state. 

As a general tendency, the effect of frustration is more pronounced in the projection 
method and it seems to be overestimated. This is already known for the static spin- 



spin correlation functions treated in the spherical symmetric approach |2l], 28 which is 
incorporated in the present calculation. This approach leads to a critical value J2/J1 = 
0.11 for the second order transition between long range order and the spin- liquid state of 
the frustrated Heisenberg model |28| which is rather low in comparison with other theories 



[H HI- 



The reasonable coincidence of the Green's function method with the exact diagonal- 
ization data at J\ = 1 is understandable since for such large values of J\ the magnetic 
energy stabilizes the size of the magnetic polaron and our concept of a polaron of small 
radius is justified. Let us note that for J\ <C t\ the bandwidth tends to be proportional 
to J\ [11-15]. In principle this result may also be obtained within our projection proce- 
dure if we enlarge the polaron size by enlarging the space of operators . To clarify the 
problem, in Fig. 2c we represent the dependence of the levels on the frustration for the 
same parameters as in Fig. 2b, but for a more realistic value of the exchange interaction 
J\ = 0.4. Comparing Figs. 2b and 2c we observe that the bandwidth for J\ = 0.4 is scaled 
from its value for J\ = 1 by a factor of 0.4 with surprising accuracy. The scaling of the 
spectrum as a whole is not as good, especially the three levels (7r/2,0), (tt,0) and (it, tt/2) 
are shifted with respect to the bottom of the band. But the order of the levels coincides 
nearly and the effects which have been discussed above may also be observed in Fig. 2c. 
Especially, the appearance of a flat band region, i.e. the decrease of the energy difference 
e(7r,0) — e("7r/2,0) is visible. So, we can conclude that the results for J\ = 1 are already 
responsible for the realistic value J\ = 0.4 if we scale the energies by a factor of roughly 
0.4. Such a scaling is already known for the case without additional hopping parameters 



7 



(*2 = *3 = 0) [H-15, 20] where the scaling is fulfilled with a slightly higher accuracy than 
may be observed here. Furthermore, comparing Figs. 2b and 2c with former results [12, 2C] 



we see that the additional hopping parameters do not increase the bandwidth. That is in 
contradiction to the work [^] where it was conjectured that the inclusion of the t2- and 
t3-terms leads to an increase of the bandwidth up to a factor of four. 

Finally, in Fig. 3 we represent the results for finite temperature T = 1 and zero 
frustration. As can be expected, in some sense the inclusion of temperature simulates the 
effect of doping by decreasing the antiferromagnetic correlation length in a similar way as 
does the frustration. But now, the exact diagonalization method is not applicable. The 
dispersion in Fig. 3 is similar to that one of Fig. lb. In particular, the extended saddle 
point between (tt/2,0) and (vr, 0) can be observed. There is one important difference that 
the saddle point in Fig. 3 lies much higher with respect to the band bottom than in Fig. 
lb. That may have consequences for the position of the Fermi level if we fill the band 
with holes. That question, however, is beyond the scope of the present investigation. 



4 Conclusion 



The ARPES experiments show different dispersions in the insulating and optimally doped 
cuprates. One observes an isotropic band minimum but no flat band in the insulating 
compound and a flat band region in the optimally doped ones. We have shown here 
that the isotropic band minimum of the undoped case may be explained due to additional 
hopping parameters in a t- J-like Hamiltonian which naturally appear in a proper reduction 
scheme from the three band Hubbard model. The doping changes the spin background 
and in the present study we have assumed that it may be simulated by a frustration term 
in the exchange energy. We have seen that frustration (doping) leads to a non-rigid change 
of the bands such that a flat band region appears. 

Let us emphasize the importance of taking into account the t<i- and i3-hopping-terms. 
As mentioned above, these hoppings lead to the isotropization of the band bottom in the 
non-frustrated case. But these terms are also important for an adequate description of the 
spectrum in the vicinity of the point (tt, 0). As it was shown in (2(J the frustration leads 
to a flat band region with an extended saddle point near the point (tt, 0) also in the case 
of the t-J\-Ji model. But then, in contrast to Fig. lb, this saddle point is extended in 
the direction (tt,0)-(tt,tt/2) and not in the direction (7r/2, 0)-(tt, 0) as in the experiment. 
So, only the simultaneous consideration of t2~, t3-hoppings and of frustration leads to the 
extended saddle point in the direction (tt/2, 0)-(ir, 0). 

In the present study we proposed to explain the different dispersion in insulating and 
doped cuprates by one characteristic set of hopping-parameters and a non-rigid band be- 
haviour. Since up to now the experiments in the insulating and doped cases have been 
done at different compounds, another theoretical approach is possible as well: different 
compounds lead to different hopping parameters and to different one-hole dispersions 
which behave more or less rigid with respect to doping. Such a point of view can be 



found in [^] and [36]. In the latter work it was proposed that the microscopic hopping 
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parameters are influenced especially by the position of the apex-oxygen. Which of the two 
approaches is the correct one can be decided in the end only by an experiment which shows 
explicitely the difference in the dispersion of the undoped and doped case of one substance. 
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APPENDIX 



The procedure to calculate the matrix elements S ab (k) ( |i~0|) and H ab (k) ( |i~2^ ) is close to 
that one given in [2C], and we present only the final expression. Let us introduce the fol- 
lowing notations for static correlation functions of Hubbard operators with noncoinciding 
indices: 

Z a = ^2 ( X i SX i+a) ) 
s 

D \ / vffSi v s±S2 yS2<j\ 

L>a,b - 2^ \ A « A i+a A i+6/ > 

i+c I 



S1S2 

T/ _ \ ^ / V~°" s l V s l s 2 V-*2«3 V s 30-\ 

Va,b,c - 2^ \ A i A «+a A i+6 A i+c / 



These expressions can be further reduced to static spin-spin correlation functions (for 
details see |2(J). A straightforward calculation gives for the overlap matrix: 

S aa = l/2 (a = 0,...,4) 

S 0a = Z a (a = l,...,4) (15) 
S ab = D a>a - b (a,6 = l,...,4) (a^b), 

and we see that it doesn't depend explicitely on the momentum k. Analogously, the 
Hamilton matrix may be expressed as: 

H{k) = K{k)+E 

K 00 (k) = 2~>^e iM , 
I 

1 i 

K aa (k) = ta Z a ^ ka + e- ika )+^ai^-a,itiV aM e m , 

l 

K ab (k) = 5^ b t b Z a e ikb + t a Z b e~ lka + t b _ a Z b _ a e ik ^ 

,a+l,a+l—b 6 ; 

I 

E = (Ji + J 2 )S + E, 
p 00 — _iv T 7 

1 m 

pOa _ _ \ " t | ^am . C n \ 

- c ' — / t "m 1^2 m ' m— a I ' 
m \ / 

paa _ _ t f ^ am 7 i r ( 7 _ ~^ m ~ 



m 

— ^ara ^fem Va—rn,a,a—b ^am $m—a,—b(Vm,a,a—b V—m,a—rn,a—m- 
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where 5 am = 1 for a = m and zero else, and S am = l — 5 am . The indices a and b (a / b) run 
over 1 to 4 and the above summation on / and m are restricted to that neighbours which 
are defined in the kinetic @ and exchange (||) parts of the Hamiltonian, respectively. 



11 



References 

[I] Tobin, J.G., Olson, C.G., Gu, C, Liu, J.Z., Solal, F.R., Fluss, M.J., Howell, R.H., 
O'Brien, J.C., Radousky, H.B., Sterne, P.A.: Phys. Rev. B 45, 5563 (1992) 

[2] Gofron, K., Campuzano, J.C., Ding, H., Gu, C, Liu, R., Dabrowski, B., Veal, B.W., 
Cramer, W., Jennings, G.: J. Phys. Chem. Solids 54, 1193 (1993) 

[3] Abrikosov, A.A., Campuzano, J.C., Gofron, K.: Physica C 214, 73 (1993) 

[4] Dessau, D.S., Shen, Z.-X., King, D.M., Marshall, D.S., Lombardo, L.W., Dickinson, 
PH., Loeser, A.G., CiCarlo, J., Park, C.-H., Kapitulnik, A., Spicer, W.E.: Phys. 
Rev. Lett. 71, 2781 (1993) 

[5] King, D.M., Shen, Z.-X., Dessau, D.S., Marshall, D.S., Park, C.H., Spicer, W.E., 
Peng, J.L., Li, Z.Y., Greene, R.L.: Phys. Rev. Lett. 73, 3298 (1994) 

[6] Friedel, J.: J. Physique 48, 1787 (1987); 49, 1435 (1987) 

[7] Tsuei, C.C., Newns, D.M., Chi, C.C., Pattnaik, PC: Phys. Rev. Lett. 65, 2724 
(1990) 

[8] Markiewicz, R.S.: J. Phys. Cond. Matter: 2, 6223 (1990) 

[9] Newns, D.M., Pattnaik, P.C., Tsuei, C.C.: Phys. Rev. B 43, 3075 (1991) 

[10] Wells, B.O., Shen, Z.-X., Matsuura, A., King, D.M., Kastner, M.A., Greven, M., 
Birgeneau, R.J.: Phys. Rev. Lett. 74, 964 (1995) 

[II] von Szczepanski, K.J., Horsch, P., Stephan, W., Ziegler, M.: Phys. Rev. B 41, 2017 
(1990) 

[12] Dagotto, E., Joynt, R., Moreo, A., Bacci, S., Gagliano, E.: Phys. Rev. B 41, 9049 
(1990) 

[13] Kane, C.L., Lee P.A., Read, N.: Phys. Rev. B 39, 6880 (1989) 

[14] Martinez, G., Horsch, P.: Phys. Rev. B 44, 317 (1991) 

[15] Eder, R., Becker, K.: Z. Phys. B - Condensed Matter 78, 219 (1990) 

[16] Hayn, R., Richard, J.-L., Yushankhai, V.Yu.: Solid State Commun. 93, 127 (1995) 

[17] Dagotto, E., Nazarenko, A., Boninsegni, M.: Phys. Rev. Lett. 73, 728 (1994) 

[18] Eder, R., Ohta, Y., Shimozato, T.: Phys. Rev. B 50, 3350 (1994) 

[19] Nazarenko, A., Vos, K.J.E., Haas, S., Dagotto, E., Gooding, R.J.: Phys. Rev. B 51, 
8676 (1995) 

[20] Hayn, R., Barabanov, A.F., Schulenburg, J., Richter, J.: One-hole motion in the 
two-dimensional frustrated t-J model, Phys. Rev. B , in press (1996) 



12 



[21] Barabanov, A.F., Beresovsky, V.M.: J. Phys. Soc. Jpn. 63, 3974 (1993); Phys. Lett. 
A 186, 175 (1994) 

[22] Jefferson, J.H., Eskes, H., Feiner, L.F.: Phys. Rev. B 45, 7959 (1992) 

[23] Belinicher, V.I., Chernyshev, A.L., Shubin, V.A.: Phys. Rev. B 53, 335 (1996) 

[24] Hayn, R., Yushankhai, V.Yu., Oudovenko,V.: Proper reduction scheme to an ex- 
tended t-J model and the hole dispersion in S^CuO^Ckj, unpublished (1996) 

[25] Hybertsen, M.S., Stechel, E.B., Schluter, M., Jennison, D.R.: Phys. Rev. B 41, 11068 
(1990) 

[26] Inui, M., Doniach, S., Gabay, M.: Phys. Rev. B 38, 6631 (1988) 
[27] Mori, H.: Prog. Theor. Phys. 33, 423 (1965) 

[28] Barabanov, A.F., Beresovsky, V.M.: Zh. Eksp. Teor. Fiz. 106, 1156 (1994) (JETP 
79, 627 (1994)) 

[29] Mermin, N.D., Wagner, H.: Phys. Rev. Lett. 17, 1133 (1966) 

[30] Novikov, D.L., Freeman, A.J., Jorgensen, J.D.: Phys. Rev. B 51, 6675 (1995) 

[31] Zubarev, D.N.: Sov. Phys. Uspekhi 3, 320 (1960) 

[32] Bacci, S., Gagliano E., Nori, F.: Int. J. Mod. Phys. B 5, 325 (1991) 

[33] Moreo, A., Dagotto, E., Jolicoeur, T., Riera, J.: Phys. Rev. B 42, 6283 (1990) 

[34] Schulz, H.J., Ziman, T.A.L.: Europhysics Letters 18, 355 (1992) 

[35] Richter, J.: Phys. Rev. B 47, 5794 (1993); Richter, J., Ivanov, N.B., Retzlaff, K.: J. 
Mag. Mag. Mat. 140-144, 1609 (1995) 

[36] Raimondi, R., Jefferson, J.H., Feiner, L.F.: Phys. Rev. B 53, 8774 (1996) 



13 



FIGURE CAPTIONS 



FIG. 1: Quasiparticle dispersion along the line M-G-X-M and contour plot for 
J 1= ti = l,t 2 = -0.08, t 3 = 0.15, T = 0.2 without frustration (a) and 
for J2 = 0.4 (b). The points in fc-space mean: M = (ir,ir),G = (0,0) 
and X = (7T, 0). The circles are the result of the exact diagonalization 
of a 4 * 4 lattice at T = with total spin Stot = 1/2. 



FIG. 2: Energy difference between the lowest state of one hole with total spin 
Stot = 1/2 for several momenta and the groundstate without any hole in 
dependence on the frustration within the Green's function method (GF) 
for Ji = 1 (a), and for the 4*4 lattice for J± = 1 (b) and J\ = 0.4 (c). 
The hopping parameters are the same as in Fig. 1. 



FIG. 3: Dispersion and contour plot for J2/J1 = and finite temperature T = 1. 
The other parameters are as in Fig. 1. 
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